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ABSTRACT 



Thermal patterns with the appearance of cyclonic motion 
and sharp thermal fronts frequently are seen in satellite 
IR images off the California coast near Pt. Sur. These 
thermal patterns are associated with distinctly structured 
upwelling systems. In nine seasonal cruises since December 
1978, nutrient fronts are strongly correlated with thermal 
fronts. However, an upwelling which is well defined by the 
satellite detected thermal front, may not be nutrient rich. 

A strong inverse linear correlation between nutrients and 
temperature though does exist for upwelling within an early 
stage of development. It is feasible, by utilizing limited 
in situ data, to infer nutrient distribution within satellite 
detected, surface thermal patterns associated with upwelling. 
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I. INTRODUCTION 



This thesis is a study of the linear correlation of sat- 
ellite detected sea surface thermal patterns to nitrate and 
phosphate levels in a coastal upwelling region. In part, it 
considers whether the chemical mesoscale in upwelling systems 
may be inferred from satellite IR detection of surface thermal 
patterns . 

The Advanced Very High Resolution Radiometer (AVHRR) on 
the TIROS (Television Infrared Observation Satellites) sat- 
ellite series can be used to locate sea surface temperature 
patterns which identify upwelling regions for ^ situ auto- 
mated biochemical and thermal sampling and analyses. Surface 
temperature and nutrient maps can be generated using ^ situ 
measurements and inference re temperature based on synoptic 
satellite IR imagery and nutrient to temperature correlations. 

Previous studies which used the same methods to detect and 
describe the thermal features (e.g., Traganza, Nestor, and 
McDonald, 1980) showed that some of these features resulted 
from upwelling pulses and produced nutrient structures of the 
same scale. These highly structured upwelling systems had 
sharp thermal and chemical gradients or fronts. The recurrent 
formation of a ’’cyclonic upwelling system” off Pt. Sur, Cali- 
fornia [Traganza, Conrad, and Breaker, 1980], offers a unique 
opportunity to investigate the significance of frontal systems in 
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regions of upwelling. Significantly higher nutrient concen- 
trations observed in upwelling systems can enhance biological 
activity in the proximity of the fronts. 

The ability to infer the nutrient concentrations associ- 
ated with an upwelling system, identified in satellite IR 
imagery, might help in explaining the relationships between 
the chemical mesoscale and pelagic ecosystem in the region. 
The occurrence of microplanktonic blooms adjacent to sharp 
thermal and chemical gradients has been documented in many 
coastal upwelling regions around the world [Barber, 1980; 
Blasco, 1980; Boyd and Smith, 1980; Traganza, Conrad, and 
Breaker, 1980]. Because phytoplankton are the primary pro- 
ducers in the pelagic ecosystem, their productivity supports 
the foodweb of which man harvests primarily the higher level 
heterotrophs . The magnitudes of nutrient concentrations and 
nutrient fluxes in upwelling systems could provide clues to 
the biological productivity in such areas. The associated 
thermo-chemical gradients or fronts could become the foci 
for fishing efforts. With the potential for biological pro- 
duction known and the growth rate predicted for an ecosystem, 
fishing fleets could be effectively dispatched to optimize 
yields . 

While thermal patterns may correlate well with satellite 
observations, correlations for chemical patterns with IR 
imagery are by inference and may be less simple, tending to 
vary as a function of the initial conditions, boundary 
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conditions, biological processes, and the dynamic processes 
acting over time. An understanding of all these processes 
may be necessary to translate satellite thermal patterns 
into nutrient levels. 

The objectives of this thesis are: 

1. To investigate the linear relationship between in 

s itu values of nutrients (nitrate and phosphate) and temper- 
ature associated with coastal upwelling occurring under 
varying nutrient and temperature characteristics of the 
ambient water mass and stages of development. 

2. To map observed values of temperature, nitrate, 
phosphate, and nutrient ratio (nitrate/phosphate). 

3. To compare the magnitudes and gradients of tempera- 
ture and nutrients associated with the oceanic front in 
different types of upwelling in different stages of develop- 
ment. 

4. To examine the relationship between the surface 
temperatures and nutrient concentrations in recently upwelled 
water and to estimate the depth from which it advected. 

5. To test the accuracy of a satellite-aided prediction 
in which the source water temperature and nutrient character- 
istics are used in a regression equation to hindcast surface 
nutrient concentrations. 
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II. METHODS 



A. STUDY AREA 

The study area (Fig. 1) extends ca. 30 km north of Pt. 

Sur and £a. 60 km to the south. It is bounded by the 122® 

30' W meridian and the Central California coastline stretch- 
ing from Monterey to Cape San Martin. 

B. CRUISES 

Three cruises were conducted aboard the Naval Postgraduate 
School's research vessel R/V ACANIA, 27 and 28 September 1979; 
29 and 30 November 1970; and 10 and 11 June 1980. Positions 
were determined by LORAN-C which has good coverage in this 
region. Running plots of temperature and nitrate and phosphate 
concentrations were maintained and used to ensure coverage of 
features of interest. Particular attention was given to the 
location of thermal gradients observed in satellite images 
received prior to each cruise. 

C. SURFACE SEAWATER SAMPLING SYSTEM 

Seawater was continuously sampled from a depth of 2.5m. 

The September and November cruises pumped seawater from a 
keel intake via the ship's pump to the shipboard laboratory. 

The June cruise used an overboard nylon hose intake at 2.5m, 
which was attached to a deck-mounted Wilden, M-4, air oper- 
ated diaphragm pump (Wilden Pump and Engineering Company, 
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Fig. 1. Bathymetry of the study area with grid superimposed. 
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Colton, California 92324) . An air debubbler was used in 
line between the pump and measuring equipment [Technicon 
Autoanalyzer, AA-II and Turner III Fluorometer; see Bronsink, 
1980] . 

D. VERTICAL PROFILE SAMPLING 

Seven Nansen casts were conducted also on the November 
cruise. Seawater samples were collected at 25, 50, 100, 200, 
300, and 400 meters for analysis of temperature, salinity, 
and nitrate. 

E. TEMPERATURE 

A thermistor sensed the temperature of the seawater 
entering the keel intake. A strip chart recorder continu- 
ously plotted it. The equipment was calibrated and monitored 
every 30 minutes by bucket thermometer readings of the sea 
surface. The mean difference, T(bucket) - T (thermistor) , of 
0.1°C (std. dev. ±0.2°C) was determined for the temperature. 
Sippican expendable bathythermograph (XBT) probes were re- 
leased every 20 minutes, or ca. every 5 km, during the Sep- 
tember cruise. The September 1979 and June 1980 cruises 
dropped XBT probes only in areas of major interest. A spac- 
ing of ca. 5 km was maintained along the major axis of the 
feature and during front transits. For calibration, sea 
surface temperature was measured by bucket thermometer 
simultaneously with the launch of each XBT probe. The mean 
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difference, TCbucket) -TCprobe), between surface measure- 
ments was 0.0*’C Cstd. dev. ±0.2°C). 

F. NUTRIENTS 

Surface (2. 5m) nutrient concentrations were analyzed 
every two minutes according to the Technicon Industrial Method 
175-72-WM, 177-72-WM [Anonymous, 1973] and 100-70-WM [Anon- 
ymous, 1978]. Nitrate here may include traces of ^ situ 
nitrite, since the nitrate is reduced to nitrite before meas- 
urement. However, according to Paulson [1972] there is little 
or no nitrite at the surface sampling depth in this area. 

G. LINEAR REGRESSION AND CORRELATION COMPUTATIONS 
Regression lines were generated for nitrate versus 

phosphate, nitrate versus temperature, phosphate versus tem- 
perature, and nutrient ratio (nitrate/phosphate) versus tem- 
perature. The linear regression analysis used the NORLSQ 
library subroutine from the W.R. Church Computer Center at 
the Naval Postgraduate School. 

The population correlation coefficient was derived by 
the equation: 

nTx^y^ - (Ex^) (Ey^) 

([nEx/- (Ex.)^[nEy/- (Ey.)M)2 

When a "none detected" condition occurred due to concentra- 
tions below the sensitivity of the instrument, the nutrient 
and temperature data were not used in the calculations. 
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H. SATELLITE IMAGERY 



The AVHRR has a temperature resolution of .5°C within 
the IR (10.5 ym to 12.5 ym) spectral window and a spatial 
resolution of 1.1 km. The thermal energy measured by the 
AVHRR is primarily a function of radiation emitted from the 
upper millimeters of the sea surface, clouds, and atmospheric 
gases. The greatest errors in estimating sea surface temper- 
ature are induced by atmospheric attenuation of the thermal 
energy due to absorption and scattering. The effect of 
absorption by moisture in the atmosphere presents the most 
severe problem because of the magnitude of absorption and the 
variability of the moisture field. At mid- latitudes over the 
ocean, this attenuation may result in measurements that are 
2 to 6°C lower than actual sea surface skin temperatures 
[Maul and Sidran, 1973]. 

The AVHRR has a temperature response from ca. -90 to 60 °C. 
The measured thermal energy is normally displayed as gray 
tones on photographic film to produce images interpretable 
as clouds, land, and water. The colder temperatures within 
these images are assigned lighter shades of gray. Image 
enhancement reduced the range of temperatures to those rep- 
resentative of the sea surface, while retaining all the gray 
tones [see Methods: Satellite Assisted Sampling Strategy]. 

In the enhancement of the TIROS-N images (Plates 1, 2, and 
3) a temperature range representative of the iui situ range 
of values was assigned to shades of gray with white equal 
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to ca. 7®C and black ca. 22®C. Studying satellite images 
alone is subjective because of the inability to assign tem- 
perature values representative of the sea surface to each 
gray shade. However, the relative temperature values and 
thermal gradients are more representative of situ measure- 
ments within the same image than are absolute values. Because 
of the atmospheric effects on the thermal energy reaching the 
satellite, these measurements can best be described as 
"apparent temperatures." 

I. SATELLITE ASSISTED SAMPLING STRATEGY 

Approximately one week prior to each scheduled departure, 
Mr. Breaker, the staff oceanographer at the National Environ- 
mental Satellite Service (NESS) , in Redwood City, California, 
would be alerted. He closely monitored images of TIROS-N 
satellite series for oceanic features of interest in the 
study area. The images were computer enhanced by him to 
better delineate sea surface thermal patterns. All images 
had the non-linear distortion associated with the earth's 
curvature and rotation removed. The day prior to the cruise 
and sometimes during the cruise, Mr. Breaker sent updated 
satellite information regarding the approximate location, 
center, size, and orientation of the feature. When possible, 
the most useful and most recent images were delivered prior 
to sailing. All proved invaluable in locating features and 
planning sampling strategy. 
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J. SATELLITE ASSISTED MAPPING OF SURFACE TEMPERATURE 



Two objective temperature fields were drawn by the CONISD 
computer program; one field was formed from only ^ situ 
temperature, the other included data from satellite imagery 
as well. 

The size of the study area was small enough that the 
spherical distortion attributed to the earth's curvature was 
assumed to be negligible. Therefore, the author constructed 
a grid on the assumed planar surface (see Appendix C) with 
the head of the Sur submarine canyon as its origin (Fig. 1). 
Each ^ situ data point had cartesian coordinates assigned. 
Using this annotated data, the CONISD library subroutine from 
the W.R. Church Computer Center drew surface contours. Be- 
cause the contouring computations involved a triangulation 
process between data points nearest the limits of the grid 
and the study area boundaries, it was necessary to assign 
ocean and near coast boundary values that were realistic to 
generate a surface map which well represented the external 
limits of the surface feature. These values were inferred 
from the thermal patterns visible in satellite IR images. 

Enhanced satellite images taken during or nearest the 
period of ^ situ observation and having a cloud-free view 
of the study area were used. These images were taken to NESS 
to geographically locate the surface thermal patterns associ- 
ated with the upwelling system on a transverse mercator 
projection of the study area. The author outlined shifts in 
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the gray tones, visible to the naked eye. These lines rep- 
resented a subjectively determined change in the gray scale. 

A uniform distribution of atmospheric moisture over the small 
area of the surface feature was assumed. The change in gray 
scale could therefore be related to synoptic variations in 
apparent surface temperature. Mr. Breaker, with his exper- 
tise in using the Baush and Lomb Zoom Transfer Scope, trans- 
cribed the apparent thermal patterns onto the transverse 
mercator projection. These patterns were then transferred 
by the author onto the grid developed for the study area. 

Using a surface temperature map generated only from in 
situ measurements, a plot of the ship’s track (to know where 
and when ^ s itu measurements were taken) and the synoptic 
thermal patterns inferred from the satellite image of the 
feature, the thermal patterns within the area of concentrated 
in situ measurements were assigned temperature values equiva- 
lent to those on the surface temperature map. The satellite 
detected and ^ situ thermal patterns matched well, within 
±.5®C. Where these thermal patterns approached the coast, 
those values inferred above were applied to the near coast 
boundary. For objective (machine) contouring the oceanic 
temperature field beyond 6 km from the seaward outline of the 

surface thermal feature was arbitrarily assigned a uniform 

0 

value. The oceanic boundary formed a curve at a constant 
distance (ca. 6 km) seaward of the surface thermal feature. 
The warmest iji situ measurement made in this oceanic region 
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vrfas assigned to the oceanic boundary. With the inferred 
near coast and ocean boundary values and ^ situ data, 
another surface temperature map was generated by the CONISD 
subroutine. Hand smoothing was applied to reduce the dis- 
tortion thought to be caused by advection or other dynamic 
processes over the time period of the study. These temper- 
ature surface maps appear in Figures 10, 20, and 34. 

K. SURFACE NUTRIENT MAPS 

The sea surface nutrient maps (Figs. 11, 12, 13, 21, 34, 
35, and 36) were generated in similar fashion to the sea sur- 
face temperature map. Each iri situ data point was assigned 
cartesian coordinates on the study area grid. The near coast 
and ocean boundary values for temperature were put into the 
regression equations for nitrate/ temperature , phosphate/ 
temperature, and nutrient ratio/temperature. The respective 
nitrate, phosphate, and nutrient ratio output were assigned 
to the coordinates of the input temperature. Using this 
annotated ^ situ and temperature inferred data, the CONISD 
library subroutine drew surface contours. Hand smoothing 
was applied to reduce the distortion thought to be caused by 
advection or other dynamic processes over the time period of 
the study. 

L. DEFINITION OF FRONT AND CALCULATION OF ITS GRADIENTS 

For the purpose of this thesis, the oceanic front is 

defined as the area where the thermal and chemical gradients 



24 



reach their maximum change per kilometer [see Discussion; 

The Relationship Between Source Water Characteristics and 
the Oceanic Front] . In order to locate the oceanic front 
and assess its gradients, the following procedures were 
carried out. 

The incremental change per kilometer of nutrient concen- 
trations and temperature was determined from the m situ 
values observed along the ship's track. The elapsed distance 
midway between the magnitudes used to compute this spatial 
gradient was assigned as the position on the ship's track at 
which the change occurred (Figs. 14, 22, and 38). This tech- 
nique was used to take advantage of the high spatial resolu- 
tion (better than 0.6 km) along the ship's track which showed 
microstructure fronts within the feature as well as the 
oceanic front. 

The magnitude of the gradient of the front was computed 
after correcting for the angle at which the ship transited 
the front. It was assumed that the front parallels isolines 
of temperature and nutriend concentrations on the machine- 
contoured charts. The angle of incidence the ship's track 
intersected these isolines at the vicinity of the front was 
measured using surface nutrient and temperature maps. Apply- 
ing the law of sines, the distance along the normal was found 
and used to compute the magnitude of the gradient. 
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III. RESULTS 



A. SEPTEMBER 1979 CRUISE 

This cruise investigated and verified the existence of a 
satellite sensed wedge- like surface thermal feature extending 
100 km SW from Pt. Sur (Plate 1). From satellite images re- 
ceived from 18 to 26 September, it appeared that this feature 
formed off Pt. Sur on 18 September. By 20 September the apex 
of this wedge- like feature had acquired a northward, anti- 
cyclonic curvature. Between 20 and 26 September this feature 
became diffuse , a new coastal upwelling event commenced 
and a new feature was superimposed. A satellite imagery 
summary is presented in Figure 2. The best image of the 
feature (Plate 1) was taken at 2255 GMT, 26 September, _ca- 2 
hours prior to the ship’s departure for the study area. 

Winds, as measured by the ship's anemometer, blew from the 
northwest averaging between 5 and 7 m* sec 

The objectives of this cruise were to obtain surface 
samples from the feature and its oceanic surroundings, to 
sample both day and night to minimize diel effects, and to 
maximize time near the cold core of the feature. A seven- 
pointed star track was planned (Fig. 3) using thermal pattern 
information derived from satellite IR imagery. Temperature, 
nitrate, and phosphate are plotted over elapsed distance in 
Figure 4. Correlation coefficients of r = 0.78 for nitrate 
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Plate 1. TIROS-N satellite IR image of the California coast 
for 26 September 1980. 
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Fig. 2. Satellite features observed from 20 September to 
2 October 1979 [Johnson, 1980] . 
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Fig. 3. Track of the September 1979 cruise and outline 
(dashed line) of the oceanic front. 
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Fig. 4. 
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Fig. 4 (cont'd) Nitrate, phosphate, and sea surface 
temperature versus elapsed distance along 
the track of the September 1979 cruise. 
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to phosphate, r = -0.91 for nitrate to temperature, and 
r = -0.66 for phosphate to temperature were obtained (Table I). 
Nutrient ratio and temperature are plotted over elapsed dis- 
tance in Figure 5. The linear correlation coefficient between 
nutrient ratio and temperature is r = -0.50. 

The linear regression analysis of nitrate/phosphate, 
nitrate/temperature , phosphate/temperature , and nutrient 
ratio/temperature (Figs. 6, 7, 8, and 9) yielded slopes of 
13.85, -2.66, -1.00, and -7.72 respectively, and x-axis in- 
tercepts of 0.45 pM phosphate, 16.68®C, 16.77°C, and 16.53°C 
respectively. The x- intercept indicated that the environment 
is nitrate limited and determined the warmest temperature for 
which these regression lines could be applied. 

The major axis of the upwelling system appeared to parallel 
and be located in the vicinity of the Sur submarine canyon 
(Figs. 10, 11, 12, and 13j . The largest area of the feature 
and the part which appears to curl cyclonically is south of 
the submarine canyon. The leading edge of the curl is located 
between the Sur and Lucia submarine canyons over the shelf 
slope, beyond the 1500 meter depth contour. 

The incremental change per kilometer of temperature and 
nutrient concentrations are plotted over elapsed distance in 
Figure 14. The oceanic front appears to be located at elapsed 
distance 13, 18, 39, 67, 166, 200, 312, 361, and 414 km 
(Table II) 
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Fig. 5. Nutrient ratio and sea surface temperature 
versus elapsed distance along the track of 
the September 1979 cruise. 
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Fig. 5 (cont'd) Nutrient ratio and sea surface 
temperature versus elapsed distance along 
the track of the September 1979 cruise. 
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Fig. 6. Nitrate versus phosphate for the September 1979 
cruise . 
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Fig. 7. Nitrate versus temperature for the September 1979 
cruise. 
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Fig. 8. Phosphate versus temperature for the September 1979 
cruise . 
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Fig. 9. 



Nutrient ratio versus temperature for the 
September 1979 cruise. 
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Fig. 10. Sea surface temperature map for the September 1979 
cruise (contour interval, 0.5®C). Map generated 
from iji situ data aided by IR imagery. 
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Fig. 11 Surface nitrate map for the September 1979 cruise 
(contour interval, 1 uM nitrate). Map generated 
from ^ situ data aided by inferences from IR 
imagery. . 
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Surface phosphate map for the September 1979 cruise 
(contour interval, 0.2 pM phosphate). Map gener- 
ated by jji situ data aided by inferences from IR 
imagery. 
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Fig. 13 Surface nutrient ratio map for the September 1979 
cruise (contour ratio interval, 1). Map generated 
by iji situ data aided by inferences from IR imagery. 
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Fig. 14 Incremental change per kilometer of temperature, 
nitrate, and phosphate versus elapsed distance 
along the track of the September 1979 cruise. 
Frontal transit occurred at peaks. 
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Fig. 14 (cont'd) Incremental change per kilometer o£ 
temperature, nitrate, and phosphate versus 
elapsed distance along the track of the September 
1979 cruise. Frontal transit occurred at peaks. 
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The thermal and chemical gradients are both sharp and 
well defined despite a narrow spread of surface nutrient and 
temperature values between the ocean and the upwelling system, 
0.01 to 4.36 yM nitrate, 0.03 to 0.96 yM phosphate, and 17.05 
to 14.80®C. The oceanic front seemed to stay in close prox- 
imity to the 15.9°C, 1.7 yM nitrate, and 0.6 yM phosphate 
isolines. The oceanic front is outlined in Figure 3. 

Figure 15 shows the vertical thermal structure of the 
upwelling system. One hundred forty XBT probes were used 
with an average spacing of 4 km. The thermal feature is 
evident from 15 to 26, 43 to 64, 197 to 246, 372 to 383, and 
408 to 428 km. The thermocline under the upwelling system 
was continuous but upwarped to within 20m of the surface. 

B. NOVEMBER 1979 CRUISE 

The staff oceanographer at NESS began monitoring satellite 
images of the study area on 9 November. He noticed a discern- 
ible though weak thermal feature on 11 November. It extended 
35 km southwest of Pt. Sur. From satellite observation, the 
thermal feature appeared to reach its maximum seaward extent 
of at least 60 km on 20 November. By 27 November it had 
become diffuse and by 28 November was indistinguishable. Due 
to intermittent cloud cover, no satellite information pertain- 
ing to a specific feature was therefore available when the 
ship departed 28 November. A prior cruise in the area from 
26 to 27 November detected a frontal feature from iji situ 
temperature measurements. 
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Fig. 15 Vertical temperature sections along the track of 
the September 1979 cruise [Johnson, 1980] . 
Vertical lines represent XBT drops. 
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Fig. 15 (cont'd) Vertical temperature sections along the 
track of the September 1979 cruise [Johnson, 1980]. 
Vertical lines represent XBT drops. 
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Fig. 15 (cont’d) Vertical temperature sections along the 
track o£ the September 1979 cruise [Johnson, 1980]. 
Vertical lines represent XBT drops. 
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15 (cont'd) 
track of 
Vertical 



Vertical temperature sections along the 
the September 1979 cruise [Johnson, 1980]. 
lines represent XBT drops. 
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As it happened, another distinctive thermal pattern was 
seen later in images on 28 November. This feature extended 
west of Pt. Sur and grew to maximum size by 6 December. It 
advected north and apparently merged with upwelling off 
Monterey Bay around 14 December (Fig. 16). Despite inter- 
mittent cloud cover, the best satellite image was taken during 
the study period at 2310 GMT, 29 November (Plate 2) . The 
average wind was 3 to 5 m*sec ^ from the northwest. 

The cruise plan had the ship transit south paralleling 
the coast. When the front was penetrated, an expanding square 
search would be followed to delineate the shape, center, 
orientation, and size of the feature. A ladder search pattern 
was then to be executed to provide additional transits across 
the major axis of* the feature. When the expanding square 
search located a small feature south of Pt. Sur, the ship 
transited north along a front and then began a ladder search 
seaward. The first leg transected the major axis of a large 
feature located west of Pt. Sur. After the general location 
of a sharp thermal gradient was found about the feature, an 
expanding square search was begun in the vicinity of the 
leading edge (Fig. 17). 

The data herein, presented is associated with those meas- 
urements taken during this ladder search and final expanding 
square. In addition, seven oceanographic stations were 
occupied near the end of the study period. Nansen casts were 
made to sample temperature, salinity, and nitrate. 
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Fig. 16 Satellite feature observed from 19 November to 
14 December 1979 [Johnson, 1980] . 
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Plate 2. TIROS-N satellite image of the California 
coast for 29 November 1979. 
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Fig. 17 Track of the November 1979 cruise and outline 
(dashed line) of the oceanic front. 
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Because of contamination in the combined working reagent, 
no reliable phosphate concentrations could be measured on 
the Autoanalyzer. 

The inverse correlation between nitrate and temperature 
CFig. 18 and 19) was strong, r = -0.93. The regression 
line of nitrate to temperature had a slope of -3.24 and an 
x-intercept of 14.74°C. The feature extended seaward in a 
southwesterly direction from a point near the coast ca- 20 km 
north of Pt. Sur. At a position ca* 20 km west of Pt. Sur, 
it appeared to curl cyclonically (Figs. 20 and 21). Unlike 
the September feature, this upwelling system did not appear 
to extend beyond the shelf slope to depths greater than 1200m. 

The incremental change per kilometer of temperature and 
nitrate are plotted over elapsed distance in Figure 22. The 
oceanic front appears to be transited at elapsed distances 
151, 196, 203, 231, 257, and 270 km. 

The thermal and nitrate gradients are sharp and well 
defined. The steepest appear to occur near the equatorward 
edge of the curl. Compared with the September cruise, the 
nitrate gradient at the front is greater for the December 
feature, whereas the temperature gradients show no great 
variation (Table II) . The spread of the surface temperature 
and nitrate values between the ocean and the upwelling fea- 
ture, 15.50 to 11.85°C and 0.12 to 9.92 pM nitrate, was 
greater in magnitude than the September cruise. The oceanic 
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Fig. 18 



Nitrate and sea surface temperature versus elapsed 
distance along the track of the November 1979 cruise. 
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Fig. 19 Nitrate versus temperature for the November 1979 
cruise . 
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Fig. 20 Sea surface temperature map for the November 1979 
cruise (contour interval 0.5®C}. Map generated 
from ijn situ data aided by IR imagery. 
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Fig. 21 Surface nitrate map for the November 1979 cruise 
(contour interval, 1 pM nitrate). Map generated 
from ^ situ data aided by inferences from IR 
imagery. 
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22 Incremental change per kilometer of temperature and 
nitrate versus elapsed distance along the track of 
the November 1979 cruise. Frontal transit occurred 
at peaks. 
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front seemed to closely approximate the 14.0®C and 2.5 pM 
nitrate isolines. The oceanic front is noted in Figure 17. 

Figure 23 shows the vertical profile of the upwelled 
feature based on 40 XBT profiles. The vertical thermal 
structure is apparent from elapsed distance 179 to 216, 226 
to 247, 275 to 290, and 298 to 319 km. The thermocline under 
the upwelling system is continuous and shows no great pertur- 
bations . 

The data obtained from the seven Nansen casts appears 
in Figures 24 and 25. Because of the increase in nitrite as 
depth increases, the actual concentration of nitrate below 
the surface layer can be less than the values recorded Csee 
Methods: Nutrients). 

Station 7 was taken near the major axis of the upwelling 
system, whereas station 2 sampled the warmer oceanic water 
to the south of the feature. Stations 3, 4, 5, and 6 were 
all taken in close proximity to the oceanic front. An in- 
trusion of warmer oceanic water into the central region of 
the feature appeared to have been recorded at station 1. The 
Nansen cast data showed that the upwelling system apparently 
corresponded . to a 26.4 tongue of water upwarping close to 
the coast and travelled seaward near the surface. The nitrate 
concentrations seem to closely agree with the density struc- 
ture. 
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Fig. 23 Vertical temperature sections along the track of 

the November 1979 cruise [Johnson, 1980] . Vertical 
lines represent XBT drops; heavy lines represent 
bottom topography. 
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Fig. 23 (cont'd) Vertical temperature sections along the 
track of the November 1979 cruise [Johnson, 1980]. 
Vertical lines represent XBT drops; heavy lines 
represent bottom topography. 
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Fig. 23 (cont'd) Vertical temperature sections along 

the track of the November 1979 cruise [Johnson, 
1980]. Vertical lines represent XBT drops. 
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Fig. 24 Salinity, temperature, and nitrate versus depth 
at seven stations on the November 1979 cruise. 
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Fig. 25 T - S and T - N diagrams of seven stations on the 
November 1979 cruise. 
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C. JUNE 1980 CRUISE 

Satellite monitoring of the waters off Pt . Sur on 9 June 
indicated a pulse of upwelled water had surfaced just before 
the ship departed Monterey en route to the study area. 

Plate 3 was the best image of the surface thermal feature, 
taken at 2242 GMT 9 June 1980, four hours prior to the com- 
mencement of the study. The wind during the study period, 

June 10 and 11, averaged between 11 and 13 m*sec ^ and was 
blowing from the northwest. 

With the feature located by satellite thermal patterns, 
a five-pointed star track was planned so that crossings of 
the feature and coastal upwelling could be made. By halving 
the ship's speed during front transits, finer resolution data 
was to be obtained to compare the poleward and equatorward 
fronts for differences in thermal and nutrient gradients. 

The cruise plan was modified because of rough seas. After 
three crossings of the feature, the front analysis was given 
priority. A modified ladder track maneuvered the ship back 
and forth across the oceanic front. Before concluding the 
study, the ship transected the upwelling along its major axis. 
The ship’s track is shown in Figure 26. 

The correlation coefficients were r = 0.96 for nitrate 
to phosphate, r = -0.96 for nitrate to temperature, r = -0.92 
for phosphate to temperature, and r = -0.96 for nutrient ratio 
to temperature (Figs. 27 and 28). Linear regression analyses 
between the above parameters yielded slopes of 18.05, -4.24, 
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TIROS-N satellite image of the California coast 
for 9 June 1980. 
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Fig. 26 Track of the June 1980 cruise and outline (dashed 
line) of the oceanic front. 
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7 Nitrate, phosphate, and sea surface temperature 
versus elapsed distance along the track of the 
June 1980 cruise. 
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Fig. 28 Nutrient ratio and sea surface temperature versus 
elapsed distance along the track of the June 1980 
cruise . 
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-0.64, and -2.07 respectively, and x-axis intercepts of 
1.07 yM phosphate, 13.47'’C, IS.SS^C, and 13.97°C respectively 
(Figs. 29, 30, 31, and 32). 

The major axis of the thermal and nutrient patterns are 
parallel to and directly over the axis of the Sur submarine 
canyon (Figs. 1, 33, 34, 35, and 36). The leading edge of 
the thermal feature appeared to curl cyclonically in the 
vicinity of the head of Sur submarine canyon. The crossing 
of the oceanic front appeared to have occurred at elapsed 
distances 14, 36, 112, 158, 172, 206, and 222 km (Fig. 37). 

The sharp and well defined thermal and chemical gradients 
at the front seem to be less than in the November feature but 
greater than in September (Table II). However, the spread of 
temperature and nutrient concentrations is the widest for the 
June cruise; 13.77 to 9.04°C, 0.29 to 19.98 yM nitrate, and 
0.91 to 2.31 yM phosphate. The oceanic front appeared to 
closely approximate the 11.5*’C, 8.8 yM nitrate, and 1.6 yM 
phosphate isolines. Figure 26 shows the oceanic front. 

The vertical temperature profile was generated from 
23 XBT profiles. The thermocline under the upwelling was 
discontinuous and appeared to surface in the vicinity of the 
oceanic front. 
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Fig. 29 Nitrate versus phosphate for the June 1980 cruise. 



72 



NITRATE. 




Fig.' 30 Nitrate versus temperature for the June 1980 cruise. 
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Fig. 31 Phosphate versus temperature for the June 1980 cruise. 
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Fig. 32 Nutrient ratio versus temperature for the June 
1980 cruise. 
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Fig. 33 Sea surface temperature map for the June 1980 

cruise (contour interval, 0.5®C). Map generated 
from in situ data aided by IR imagery. 
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Fig. 34 Surface nitrate map for the June 1980 cruise 
(contour interval, 1 pM nitrate). Map 
generated by ^ situ data aided by inferences 
from IR imagery. 
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Fig. 35 Surface phosphate map for the June 1980 cruise 
Ccontour interval, 0.2 yM phosphate). Map 
generated from ijx situ data aided by infer- 
ences from IR imagery. 



78 




Fig. 36 Surface nutrient ratio map for the June 1980 cruise 
(contour ratio interval, 1). Map generated from 
in situ data aided by IR imagery. 
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Fig. 37 Incremental change per kilometer of temperature, 
nitrate and phosphate versus elapsed distance 
along the track of the June 1980 cruise. Frontal 
transits occur at peaks. 
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LINEAR REGRESSION ANALYSIS 
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TABLE II 



GRADIENTS AND OBSERVED VALUES AT THE OCEANIC FRONT 
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IV. DISCUSSION 



A. SATELLITE DETECTED THERMAL PATTERNS COMPARED 

TO ^ SITU SURFACE MAPS 

Inferring the distribution of surface nutrient concen- 
trations from surface thermal patterns detected by satellites 
appear feasible from the correlation of ^ situ nutrients and 
temperature. In this study sea surface nutrient maps were 
generated for comparison with satellite detected thermal 
patterns formed by upwelling systems off Pt . Sur, California. 
For example, the surface nitrate and phosphate maps for the 
June 1980 feature (Figs. 34 and 35) appear similar to the 
thermal distribution (Fig. 33) with regard to the general 
location of the nutrient feature, its center, and the orien- 
tation of the major axis. However, the nutrients exhibited 
a highly structured pattern with sharp gradients within the 
feature which was not predictable from the satellite image. 
The magnitude of the difference between the distribution 
implied by the image and the actual surface distribution 
of nutrients within the feature varied among the April, Sep- 
tember, November 1979 and June 1980 upwelling systems. (The 
other features lacked sufficient data to construct surface 
nutrient maps.) The thermal and spatial resolution of the 
AVHRR (1°C and 1.1 km respectively) could account for some 
of the difference. However, the relative ranges of the 
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nutrient values compared to the range of the temperature 
varied too much in different features at different times to 
be satisfactorily explained by the lack of resolution alone. 

B. NUTRIENT-TEMPERATURE CORRELATION AND 

"ACTIVE UPWELLING" SYSTEMS 

To better evaluate these differences, regression analyses 
were applied to ^ situ values of nutrients and temperature 
within the different features which formed at different times 
of the year (Fig. 38). Although there were variations, the 
slopes for all features were negative. A strong inverse 
linear correlation existed in five of the seven features 
observed since December 1978 (Table I) . This implied that 
regression lines between nutrients and temperature might be 
useful in describing the distribution of nutrients with re- 
spect to surface temperature patterns. The variations in 
slope of these regression lines suggests that the relation- 
ship of nutrients to temperature might depend on numerous 
factors. These may include the temperature and nutrient 
characteristics of the source water, the depth from which 
the water was advected, wind stress and bathymetric effects, 
surface dynamic processes (viz., mixing, advection, and heat 
exchange) and biological processes. Because of the complexity 
of the problem, trends were investigated with a view toward 
developing indices for prediction of nutrient distribution. 

The strong inverse correlations between temperature and 
nutrients occurring in five cruises suggested that ' active 
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NUTRIENT RATIO ‘ NOj / PO4 PHOSPHATE, NITRATE. 



CRUISE OATES 




SEA SURFACE TEMPERATURE. *C 



Fig. 38 Summary of regression lines for nitrate, phosphate, 
and nutrient ratio versus temperature 7 DEC 1978 - 
10 JUN 1980. 20 JAN, oceanic water; 7 DEC, 30 APR, 

8 AUG, and 10 JUN upwelled water from below thermo- 
cline; remaining cruises "upwelled" water from 
above thermocline. 
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upwelling' systems had occurred. This was readily apparent 
for the April 1979 and June 1980 cruises for which the ver- 
tical temperature profiles showed the thermocline surfacing 
near the edge of the surface thermal feature. This break, 
through the thermocline, would have permitted the advection 
of nutrient-rich, "biochemically new" water [Traganza e^ al . , 
1980], to the surface waters within the feature. This seems 
to be supported by the high nutrient concentrations, high 
nutrient ratios, and low temperatures below 10®C observed 
inside the feature (Figs. 27 and 28, respectively). 

Although satellite history of the November 1979 cruise 
suggests the observed feature was sampled during a stage of 
'active upwelling ', XBT data indicate that the s ystem appeared 
not to have penetrated the thermocline. The thermocline under 
the feature was continuous and did not appear to be perturbed 
by the shallow (ca. 50m) upwelling system. This would have 
accounted for the moderate level of nitrate observed at the 
surface inside the feature. 

The December 1978 and June 1979 cruises did not collect 
vertical temperature data. From the wind data and the mod- 
erate levels of nitrate and phosphate, the upwelling system 
observed on the June 1979 cruise was assumed ' active ' 

[Conrad, 1980]. Because the nutrient levels were similar 
in magnitude to the November 1979 cruise, this feature could 
be representative of a shallow upwelling. 
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As the December 1978 cruise had low temperatures, high 
nutrient concentrations and nutrient ratios similar to those 
of the April 1979 and June 1980 cruises, the feature observed 
could have been associated with a discontinuity in the ther- 
mocline permitting deep upwelling. Thus, strong linear 
correlations between nutrients and temperature may be assumed 
for 'active upwelling' systems. 

C. THE RELATIONSHIP BETWEEN SOURCE WATER 

CHARACTERISTICS AND THE OCEANIC FRONT 

The ability to determine the initial condition of the 
upwelled water was considered a major factor in prediction. 

The nutrients in this subsurface source water would provide 
the maximum nutrient concentrations observed in the upwelled 
surface water. The nutrient concentrations and temperature 
were expected to be a function of depth and ambient water 
mass. The depth of the source water is most likely tied to 
the duration and strength of wind or other forcing functions, 
e.g., internal waves [Johnson, 1980]. 

Within an ' active upwelling, ' the surface is replenished 
by vertical advection of source water. The properties char- 
acteristic of the source water are confined to the surface 
feature by the oceanic front. For every feature studied, the 
maximum gradient analysis of ijn s itu temperature and nutrients 
(see Methods: Definition of Front and Calculation of its Grad- 
ients) showed the nutrient fronts coincided within +.6 km to 
the thermal front (Figs. 14, 22, and 37). Thus, the thermal 
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front detected in satellite imagery closely approximates 
the nutrient front. 

Since nutrients are linearly related to temperature, 
thermal gradients are well represented within a satellite 
image, and the nutrient fronts and thermal front coincide, 
the features were examined to see if the sharpness of the 
thermal front might reflect the range of nutrients within 
the feature. 

The sharpness of the thermal front was not found to be 
a dependable measure of the nutrient distribution within a 
feature. The thermal gradients at the fronts of the December 
1978, April and September 1979, and June 1980 features had 
similar magnitude, 0.12®C/km to 0.38°C/km (Table II). Within 
the December, April, and June features high concentrations, 
8.3 to 23.1 yM nitrate and 1.0 to 2.3 yM phosphate were ob- 
served. These surface values were considered the result of 
'active' deep upwelling. However, the September feature had 
significantly different nutrient concentrations, 0.6 to 4.3 
yM nitrate and 0.5 to 0.9 yM phosphate. The vertical temper- 
ature profile of the September cruise showed a continuous 
though greatly perturbed thermocline under the surface feature 
(Fig. 15). The low nutrient concentrations and nutrient 
ratios, together with the vertical temperature date [see 
Johnson, 1980] , suggest that a shallow upwelling or mixing 
in the upper layer caused the feature. The satellite history 
of the feature revealed the thermal patterns were becoming 
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indistinguishable at the time of the cruise. The lack of 
strong correlation in the nitrate to phosphate and nutrient 
to temperature regressions also suggested that the feature 
was dissipating. The June, August, and November 1979 fea- 
tures also showed differences between nutrient distributions 
despite their similar thermal fronts. 

D. NUTRIENT DISTRIBUTION PREDICTION REQUIRES GROUND TRUTH 

Because of the inability to describe adequately and there- 
fore predict the nutrient distribution from satellite detected 
thermal patterns alone, ground truth data are needed. In 
that the linear correlation was strong between nutrients and 
temperature, an indicator of active upwelling, for the five 
features mentioned earlier, the sum effects of the initial 
conditions, boundary conditions, and dynamic processes appear 
confined to altering the slope of the nutrient to temperature 
regression line and its intercept. 

E. THE FEASIBILITY OF SATELLITE PREDICTION 

AIDED BY LIMITED ^ SITU SAMPLING 

A test, therefore, was conducted to determine how good a 
satellite prediction could be, using source water temperature 
and nutrient characteristics in a regression equation to 
forecast surface nutrient concentrations. The regression 
line is bounded by the coldest temperature in the feature 
and the warmest oceanic water temperature. 
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1 . Use of the Oceanic Water Temperature as an 
Intercept for a Predictive Regression Line 

The temperature CT) intercepts computed from ^ situ 
nutrient/temperature regression equations were compared with 
the warmest oceanic temperature to determine if a trend could 
be seen that would be useful in developing a regression equa- 
tion to forecast nutrient concentrations. The phosphate/ 
temperature T- intercept was always greater (warmer) than the 
nitrate/temperature T-intercept. Because the oceanic en- 
vironment appeared to be nitrate-limited (nitrate goes to 
zero while phosphate is still present, see Figs. 6 and 29), 
the phosphate/temperature line would have to be extrapolated 
to reach the temperature axis, thus greater (warmer) temper- 
atures. This would be an over estimate of the actual oceanic 

water temperature. The difference, V^an ' intercept > 
active upwelling systems as defined by nutrient values and 
concentrations ranged from -0.21 to +0.61°C for nitrate/ 
temperature and -0.18 to +0.27°C for phosphate/temperature 
regression equations (Table I) . Within the thermal resolu- 
tion (.5®C) of the satellite sensor, the T-intercept might 
therefore be considered equivalent to the temperature of the 
oceanic water. 

2 . Hypothesis That Nutrient/Temperature 
Relation within Source Water Is Linear 

Strong linear correlations were seen between both 

nutrients and temperature and their gradients in the November 
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1979 shallow upwelling system and June 1980 deep upwelling 
system. From satellite history, these features appear to 
have been observed during an initial stage of development 
indicating the water upwelled from its source depth had been 
in contact with the surface for a day or less. Therefore, 
the source water nutrient and temperature characteristics 
should be linearly related. If this is the case, and the 
source water regression line is known (from recent or his- 
torical data), this relation could be used, along with the 
oceanic water temperature as a T-intercept, to provide a 
first approximation of the nutrient distribution reflected 
by the surface thermal patterns. 

In order to test this hypothesis, source water char- 
acteristics were examined using T - S and T - N diagrams from 
the shallow upwelling November 1979 Nansen casts (Fig. 25). 

The 2 and 25m T - N samples for station 7 (located near the 
coldest water in the feature) and stations 3, 4, 5, and 6 
(located near the front) were plotted on the T -N diagram for 
station 2, which was representative of the oceanic environment 
(Fig. 39). The 25m samples were tightly clustered on the 
station 2 curve between the 25 and 50m depths, suggesting 
that the water had upwelled from that depth. The 2m samples 
of the ocean front stations were skewed off the station 2 
reference curve toward warmer, lower nutrient water. However, 
station 7 varied little between the surface and 25m depth. 
Apparently station 7 was nearest the area where the source 
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TEMPERATURE, ‘C 



Fig. 39 T - S and T-N diagram for station 2; ocean reference, 
with 25m (solid circles) and 2m (open circles) data 
from stations 3, 4, 5, 6, and 7 added to identify 
source water. 
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water first came in contact with the surface. The ocean 
front station data were obtained in this source water after 
it had been affected by mixing, heat exchange, and other 
processes, dynamical and biological. 

3 . The Effects of Wind Mixing and Heat Exchange to 
the Atmosphere on Surface Water near Oceanic Front 

The difference in 2 and 25m temperatures could be 

associated in part with the mixing caused by the 3 to 

5m* sec ^ winds during the November cruise and the atmospheric 

thermal exchange associated with a low percentage of cloud 

cover and £a. 10 hours of daylight observed. Over the one 

day, the upwelled water known to have been in contact with 

the surface the longest (which was assumed to be the water 

at the oceanic front) was predicted to be ca. 0.8°C warmer 

[James, 1966]. This agreed with what was observed on the 

T - N diagram (Fig. 39) . 

4 . Determination of Depth From Which 
Source Water Was Advected 

It was postulated that a nutrient source depth range 
could be determined by the depths at which the cold tempera- 
ture in the center of the feature, T^, and the temperature 
of the ocean front, T^ , were found on a temperature versus 
de.pth profile, obtained, e.g., with an XBT. From these depths 
a range of nutrient values could be measured by ^ situ 
sampling. The change in nutrients versus the change in 
temperature over the depth range, corresponding to - T^ , 
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could be a first approximation of the slope, m, for a regres- 
sion equation to predict nutrient concentrations within the 
surface feature. This slope approximation could be combined 
with the oceanic water surface temperature as described 
earlier to compute the nutrient CN) intercept, b, whereby 
the total regression equation, N = mT + b, can be formulated. 
The conditions for which this approximation would be best 
suited are those in which there were not thermal exchanges 
with the environment and changes in nutrient concentrations 
due to biological activity. Accuracy of the application to 
IR satellite data would depend on the extent of agreement 
between the surface temperature inferred from the satellite 
image and the actual surface temperature. 

\ 

5 . A Hindcast of Nitrate'* Concentration Using 

Data from the November 1979 Cruise 

The November 1979 cruise station 2 vertical tempera- 
ture profile showed the nitrate source depth range for 
11.8S°C, T^, and 14.0®C, T£ , to be 40 to 14m respectively. 

The nitrate range between these depths was 11.7 to 3.9 pM. 

The slope for this line, therefore, was c^. -3.7. A second 
approximation correcting for wind mixing and thermal exchange 
with the atmosphere [James, 1966] used 13.2®C for T^ to gen- 
erate a slope of ca. -3.1. The oceanic water surface temper- 
ature used to predict the N-intercept, b = -mT, was 15.50°C. 
From these few measurements, a predictive regression equation 
was formulated. 
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The predictive regression equation computed from 333 



in situ measurements was compared with the equations developed 
using the slope derived from source water characteristics and 
the temperature of the ocean water used as the T- intercept. 

The nitrate concentrations for 11.85® and 14.0®C were computed. 

The best fit linear regression (correlation, r=-0.93) 
based on underway iji s itu measurements was: 

N = -3.24 T + 47.74. 

where N = pM nitrate and T = ®C temperature. The predictive 
regression equation, assuming no thermal exchange with the 
environment, was: 

= -3. 7 T + 57.35 

The second approximation which corrected for wind mixing and 
atmospheric thermal exchange was: 

N 2 = -3.1 T + 48.05 

a. Comparison of Predicted Nitrate Concentration 
to That Computed by Regression Analysis of 
• Observed Data 

The solution to the equations are as follows: 

for T^ of 11.85®C, N = 9.82 pM, = 13.5 pM, and N 2 = 11.3 pM 

nitrate; for T. of 14.0°C, N = 2.94 pM, N, = 5.5 pM, and 
^ N, - N 

= 4.7 pM nitrate. The percentage of error ( x 100 

^ N 2 - N N 

and — — X 100), associated in predicting the nitrate 

concentration at the cold center or most recently upwelled 

water appears to be ca. +37% for the first approximation. 
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N^, and ca. +15% for N 2 which accounted for wind mixing and 
thermal exchange. In the vicinity of the ocean front, the 
percentage of error increased; for of 14.0°C, had an 
error of ca. +98% whereas N 2 had an error of C£. +58%. In 
both hindcasts the nutrient concentrations predicted were 
higher than those computed from the equation based on iji situ 
surface measurements. Both hindcasts also deteriorated in 
reliability proceeding away from the cold center of the fea- 
ture. The predictive equation which corrected for wind mixing 
and thermal exchange with the atmosphere had significantly 
less error than the equation that did not. Because nitrate 
and temperature are nonconservative within the upper layer 
of the ocean, the 15% error derived from linear regression 
should be considered reasonable. In all, the longer the up- 
welled water is in contact with the surface layer of the 
ocean, the more complex the interactions with the physical 
and biological environment may become, and the more unreal- 
istic the simple predictive equations. 

b. Summary of Conditions for Which This 

Test Appeared Satisfactory 

This hindcast was applied only to a shallow up- 
welling system which satellite history showed to be in an 
initial stage of development. The prediction was based on 
knowledge of the temperature of the sea surface within the 
coldest thermal pattern detected, surface temperature of the 
oceanic front detected by thermal patterns, a vertical 
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temperature and nutrient profile within the upper 200m of 
the ocean seaward of the oceanic front, and atmospheric 
parameters to hindcast ocean-atmosphere thermal exchange 
[James, 1966]. In this special case, given this limited 
ground truth, it was possible to hindcast within 151 error 
the maximum nutrient concentration in an upwelling system 
detected by satellite imagery. It therefore seems reasonable 
to assume that at least the major patterns of nutrient dis- 
tributions can, with further research, be inferred using 
satellite imagery and limited ground truth (i.e., data from 
buoys and AXBT's). 
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V. CONCLUSIONS 



1. Inference of nutrient distribution by satellite 
detected upwelling systems is feasible. 

2. Active upwelling systems are expected to have strong 
inverse linear correlations between nutrients and temperature. 

3. The nutrient front position can be approximated 
closely by the thermal oceanic front. 

4. The nutrient distribution within a feature can not 
be related to the sharpness of the thermal front. 

5. To predict nutrient distributions, ground truth as 
well as satellite detected thermal patterns are required. 

6. A linear regression can be used to forecast nutrient 
maxima for upwelling systems in the initial stage of develop- 
ment aided by only limited ^ situ data. 

7. The approximation of nutrient concentrations by linear 
regression can be improved by estimating the effects of wind 
mixing and thermal exchange with the atmosphere. 

8. With greater knowledge of source water character- 
istics (from ^ situ monitoring or historical data) , stage 
of development (inferred from satellite images and ^ situ 
monitoring), and dynamic processes (wind mixing, advection, 
and heat transfer) a forecast of nutrient distributions with 
a surface thermal feature could be made. 
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APPENDIX A 



Listing of Surface Data: Time, Latitude, Longitude, 

Elapsed Distance, Nitrate, Phosphate, Nutrient Ratio, 
ATP, AATP/ATP, Chlorophyll, Temperature 
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APPENDIX B 
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